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We present analytic self-similar or traveling wave solutions for a one-dimensional 
coupled system of continuity, compressible Euler and heat conduction equations. 
Different kind of equation of states are investigated. In certain forms of the equation 
of state one can arrive to a picture regarding the long time behavior of density and 
pressure. The impact of these quantities on the evolution of temperature is also 
discussed. 

PACS numbers: 34.10,+x, 34.50.-s, 34.50.Fa 



I. INTRODUCTION 

Studying hydrodynamical systems with heat conduction is an important task, both from 
application and from theoretical side as well. There are numerous monographs and textbooks 
available which investigate three dimensional real fluid mechanics with heat transfer from the 
engineering point of v,ew Q. UsuaUy, all the technical q nestion 8 lihe the applied nn m enc S , 
grids are analyzed in depth as well. 

In the following we will investigate the most simple one dimensional system of continuity, 
compressible and non-viscous Euler and heat conduction equations. We try to find analytic 
results, as a first point self-similar physically important diffusive solutions are investigated. 
If such solutions are not to find, as a second point traveling wave solutions are obtained and 
discussed. 

These are basically two distinct very powerful trial functions to investigate the global 
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properties of the solutions of various single non-linear partial differential equations 2|-[5]. 
One is the so-called self-similar solution which describes the intermediate asymptotics of a 
problem: they hold when the precise initial conditions are no longer important, but before 
the system has reached its final steady state. They are much simpler than the full solutions 
and so easier to understand and study in different regions of parameter space. A final reason 
for studying them is that they are solutions of a system of ordinary differential equations 
and hence do not suffer the extra inherent numerical problems of the full partial differential 
equations. In some cases self-similar solutions helps to understand diffusion-like properties 
or the existence of compact supports of the solution. 

The other method is the traveling wave Ansatz which gives us a deeper insight into the 
wave properties of the system like the propagation speed which can be even time dependent. 

In our former work the paradox of heat conduction was investigated with the self-similar 
Ansatz [6j and gave a new equation which has a finite propagation speed. Recently, a study 
about the generalized Cattaneo law is under publication where both the heat conduction co- 
efficient and relaxation time have power law temperature dependence 7j]. The results can be 
continuous or shock wave-like with compact support. With a straight forward generalization 
of the self-similar Ansatz even a partial differential equation(PDE) system the 3 dimensional 
Navier-Stokes equations was successfully investigated [8-]. Gilding and Kersner [5J studied 
large number of nonlinear diffusion-convection problems with the help of traveling waves. 

As a beginning, without completeness we mention some analytic studies about hydrody- 
namical systems which are available in the literature. 

Beyond traveling waves and self-similar Ansatze there are various examination techniques 
available in the literature. Manwai js] studied the N-dimensional (N > 1) radial Navier- 
Stokes(NS) and Euler equation with different kind of viscosity and pressure dependences 
and presented analytical blow up solutions. His works are still 1+1 dimensional (one spatial 
and one time dimension) investigations. Another well established and popular investigation 
method is based on Lie algebra by which numerous studies are available. Some of them are 



even for the three dimensional case, for more see |lfl. Unfortunately, no explicit solutions 
are shown and analyzed there. Fushchich et al. |ll( construct a complete set of (5(1, 3)- 
inequivalent Ansatze of codimension 1 for the NS system and they present 19 different 
analytical solutions for one or two space dimensions. They last solution is very closed to our 
former one js] but not identical. Further two and three dimensional studies of the Navier- 



Stokes equation based on group analytical method were presented by Grassi [12]. They also 
presented solutions which look almost the same as ours |8] , but they consider only two space 
dimensions. 



Recently, Hu et al. 13j presents a study where symmetry reductions and exact solutions 



of the (2+l)-dimensional NS were presented. Aristov and Polyanin 14j] use various methods 
like Crocco transformation, generalized separation of variables or the method of functional 
separation of variables for the NS and present large number of new classes of exact solutions. 
Sedov in his classical wo* [2] (Page 120, presented analytic solutions for the three dnnen- 
sional spherical NS equation where all three velocity components and the pressure have 
polar angle dependence (9) only. Even this kind of restricted symmetry led to a non-linear 
coupled ordinary differential equation system which has a very rich mathematical structure. 
Some similarity reduction solutions of the two dimensional incompressible NS equation was 
presented by Xia-Ju [ijj]. Additional solutions are available for the 2+1 dimensional NS also 
via symmetry reduction techniques by |16| . There is a full three dimensional Lie group anal- 
ysis is available for the three dimensional Euler equation of gas dynamics, with polytropic 
EOS unfortunately without any heat conduction mechanism [it}]. Of course one may find 
numerical methods for solving equations of fluid dynamics The topic is quite important, 
there are interesting applications of the equations of fluids in granular matter 

To our knowledge there are very few analytic (self-similar or traveling wave) solutions 
known for any non-linear partial differential equation systems (PDEs) till today [^J. Our 
experience shows that there are no such comparative studies available for fluid mechanics 
with heat conduction in the literature which we present on the next pages. 

In the next section we outline our starting model in the later subsections we give various 
solutions for three different polytropic equation of states (EOS), plus for the virial EOS. If 
our results meet some other theories (like porous media equation) then the corresponding 
literature is mentioned. The paper ends with a short summary. 



II. THE BASIC MODEL 



We start with the following one-dimensional compressible and non-viscous fluid with 
ordinary heat conduction which means the following coupled PDE system of continuity, 
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Euler and heat conduction 

p(x,t) t + [p(x,t)v(x,t)) x = 0, 

v(x, t) t + v(x, t)v(x, t) x = -^—-p(x, t) x , 

p(x,t) 

T(x, t) t + v(x, t)T(x, t) x = XT(x, t) xx , (1) 

where p, v, T, P are the density, velocity, temperature and pressure field or distribution, 
respectively. Subscript means partial derivation with respect of time and the x to the 
Cartesian coordinate. We skip viscous velocity term in Euler and in the energy equation 
as well. As a starting point we include the polytropic EOS p(x,t) = ap(x,t) n were n is a 
real number and a is a material constant and A is the heat conduction coefficient. After the 
investigation of this basic equation system we will introduce modifications for the equation 
of state. 



In the following we are__ 
which is well-known from 



ooking for the solution in the form of the self-similar Ansatz 
2r|4| 

V(x,t) = r a f(^) :=r-/fo) (2) 

where V(x, t) can be an arbitrary variable of a PDE and t means time and x means spatial 
dependence. The similarity exponents a and /3 are of primary physical importance since 
a represents the rate of decay of the magnitude V(x,t), while j3 is the rate of spread (or 
contraction if /3 < ) of the space distribution as time goes on. The validity of this self- 
similar Ansatz also means that the solution has no characteristic time scale. The most 
powerful result of this Ansatz is the fundamental or Gaussian solution of the Fourier heat 
conduction equation (or for Fick's diffusion equation) with a = (3 = 1/2. The function f{rf) 
is called the shape function. 

For our system we use the following notations and shape functions 

T(x, t) = ry , v(x, t) = r s g (^) , p(x, t) = t^h (3) 

where the new variable is rj = x/t^. Calculating all the first time and spatial derivatives of 
03]) and putting them back to ([!]) after some algebra we get the following non-linear ordinary 
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differential equation (ODE) system 

- 1 h- -r]h! + gti + hg' = 0, 
—j9 - + 99 = ~ a I — — J hK 1 )h , 

-af - \nf + gf = Xf", (4) 

where prime means derivation with respect to rj. Among the initially free parameters 
a,/3,7,n we got the following constraints a, 7 are still arbitrary but (3 = 5 = 1/2 and 
n = (l+7)/7- The material constants c, A are still independent. Note, that the first and 
second equations of ([I]) are independent of the third one. 



A. Aspects coming from the conservation law of density 

The first equation is a conservation equation (without source term) and if 7 = 1/2 it can 
be integrated 

' //,V = W (5) 



2 

giving g = rj/2 and g' = 1/2 when h ^ and the integration constant was chosen zero. 
Plugging this result back into the second equation of ([I]). 

-1 = Sahti (6) 
4 

where n = (1 + 7V7 = 3. Which means that only the p = ap 3 ploytropic EOS can give 
self-similar solutions in this flow system. The solution of ((6]) can be obtained via a simple 
integration giving the following expression 



3a 

where C\ is the integration constant (we will set it to 0) and a is still the material constant 
of the polytropic EOS. Last, let's investigate the third heat conduction equation, with the 
knowledge of the velocity its solution becomes quite simple 

-af = Xf". (8) 

The solution is well known 

far)\ . (oiT]\ , . 

/ = c 2 cos [-J- J + c 3 sm \— j (9) 
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At last, we summarize all our results 



9= 2> k 



4rj 2 + 2ci 



3a 



f = c 2 cos J + c 3 sm 



(10) 



where n = 3, ^ = (3 = 8 = 1/2 and a is arbitrary. Now the velocity, density and the 
temperature field read 



v(x,t)=r d g( V ) = - 



p(x,t) = r^hirj) = rH 



T(x,t) = t- a g(r]) = t~ a c 2 cos 
and the applied polytropic EOS 

p = ap n = a{r^h{r])) 3 = a 



ax 



fe) 2 + 2 ^ 

3c 
/ ax \ 

+C3Sm (w = t)_ 



(11) 



t"5 



3c 



(12) 



All the solutions have physical meaning for correspondingly bounded space and sufficiently 
large times. We will see later on that according to the continuity equation the velocity field 
is always a simple function of time and coordinate, however the temperature distribution 
have a much complexer form. In this case it is a quickly oscillating and decaying function 
or a slowly oscillating and slowly decaying function. 



B. Aspects on solution starting from the equation regarding the velocity 

We had the possibility to see that in certain cases, when the first equation is integrable, 
solutions of the system can be found. This was accessible when the pressure is proportional to 
the third power of the density. Now we try to find solutions when the pressure is proportional 
to the second power of the density 



v=- 2 p 



(13) 



The role of b is to fix the proper physical dimension. In this case the second equation of ([I]) 
can be written 

v t + vv x = -bp x (14) 
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If we consider that the system is not too far from the steady state, correspondingly the 
change in time of the gradient of the density is negligible, we get 

v = —bp x -t + c, (15) 

where c can be related to initial velocity and initial gradients at time to if these terms do 
not cancel each other. We suppose that ( |T5i) holds for a longer time scale than the gradient 
would undergo to considerable changes. We are interested in a density profile in such a 
quasi- stationary state. The procedure can be also regarded as taking a first initial condition 



21] . Inserting (j!5p into the mass conservation 



of the velocity field and iterate it successively 
formula we get 

Pt - tb(p x ) 2 + cp x - tbpp xx = 0. (16) 
If c = and dividing the equation with t we obtain the following equation 

^ - b( Px ) 2 - bp Pxx = 0. (17) 

Dividing the last equation by b and multiplying by 2, we have 

p tl - 2{p x f - 2p{ Pxx ) = (18) 

where t\ = bt 2 /A. This equation in certain aspects may resemble to the Kardar-Parisi- 
Zhang(KPZ) equation [22] which describes the dynamics of surface growth. The original 
KPZ equation (where the third term doesn't contain the p variable) was investigated with 



the self-similar Ansatz by the recent authors [23[ and found physically reasonable solutions 
which can be expressed with the error function. 

The last equation can be brought to a more compact form 

Pt, = (P 2 )xx, (19) 

There is an even more important relationship with the nonlinear heat conduction or porous 
media equation which has the form of 

fhx = (P m )xx, m>l. (20) 



In 



24J | Zeldovich and Kompaneets have found the fundamental solution pi of this equation 



which we write in the following form: 
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where A is constant and the subscript + at the bracket indicates that only physically relevant 
- positive - solutions are taken into account. Regarding the other constants we have 

a = f3 = —^— } B 2 = m ~ l - . (22) 
m + 1 2m(m + 1) 

One can see that this solution has bounded support in x for any t\ > which is a hyperbolic 
property. Using comparison principle for such equations one can show this finite speed 
property for any initial condition having compact support. However, the fronts are not 
straight lines: x(ti) = , /3 < 1 so the speed of propagation x(tx) goes to zero if t\ goes 
to infinity. One can also see that p\ is of source-type: pi(x, 0) = Ki8(x). 

In our case m = 2, consequently a — /3 — 1/3, and the solution for the density distribution 
reads 

tl ) = t - 1/3 (A 2 - B 2 x% 21 *) = V3 r 2 / 3 ||a 2 - 5V ^ ^ t- 4 / 3 j . (23) 

This solution gives us a density profile for a bounded spatial region and after a sufficiently 
long time. The velocity field is obtained via (15) 



v(x,t) = ^. (24) 



For the temperature distribution the 



T t + ^T x = \T XX (25) 

PDE have to be solved. Considering the T(x,t) = i~ 7 / (jH the following ODE is derived 

Xf" " ^ + 7/ = (26) 

with arbitrary 7 and for u = 1/2. The solutions are the Kummer M and Kummer U 
functions of the form of 

1 ~ 3 1 2 \ (\ 3 1 2 



/ = Cl M [ - - -, — „ + c 2 U - - 3 7 , -, — „' ) (27) 



Exhausted mathematical properties of the Kummer function can be found in 
completeness we give the temperature distribution as well, 



25j. For 



T(x,t) = r 1 



' (\ 3 x 2 \ x TT (\ 3 x 2 \ x 



(28) 
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In contrast to the first case, we applied a simplification in the Euler equation (14) therefore 
we can consider solutions which are outside the scope of the self-similar class. Therefore, we 
may try to find other physically important solutions, like traveling-waves with the Ansatz 
of p = h(() = h(x — ct) where c is usually the wave propagation velocity. (We still use h as 
shape function). However, if we consider (1171) . it can be rewritten in the following form 



Pt 2 ~ KPxf - bpp xx = 



where t 2 = t 2 /2. We take the 



p = h(() = h(x - ct 2 ) = h(x - at 2 /2) 



(29) 



(30) 



Ansatz where instead of the constant propagation speed c one may rather consider a kind 
of constant acceleration, and the constant c is renotated with a. We do this in the purpose 
of finding analytical solutions, and to get the results closer to a physical interpretation. 
Inserting this formula into (j29|) we obtain the following ODE: 



bhh" + h'{bh! + a) = 



(31) 



The solution is the following 



MO 



bc\ 



LambertW 



b 2 c 1 



bc\ 



+ 1 



(32) 



where Ci,C2 are integration constants and the Lambert- W function, (which also called the 
omega function), is the inverse function of 



f(W(x)) = W(x)e w (x) 



(33) 



and can be evaluated with the following series expansion 

W(x) 



n=l 



(n- 1)! 



(34) 



Additional properties of this function can be found in [26] . Banwell and Jayakumar [2j 
showed that a W-function describes the relation between voltage, current and resistance in 
a diode, and Packel and Yuen 28) applied the W-function to a ballistic projectile in the 
presence of air resistance. Other applications have been discovered in statistical mechanics, 
quantum chemistry, combinatorics, enzyme kinetics, the physiology of vision, the engineering 
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of thin films, hydrology, and the analysis of algorithms [29]. Note, that in contrast to other 
special functions (e.g. Whitakker, Bessel, Kummer functions) the LambertW function is not 
so widely used or applied. The final form of the density field is 



p(x - at 2 /2) = ^ 



a 



e b2 °i 

LambertW | ; | + 1 

bc\ 



(35) 



Figure 1 presents the density filed of ( 155]) for the ci=c 2 = a = b= l parameters. Note, 
that the solution is a steep wave front. 

Inserting the space derivative of the density into ( JT5|) and multiplying with the time the 
velocity field reads the following 



LambertW — 



feci 

v{x - at 1 j2) = (36) 

LambertW I ^ ) + 1 

Figure 2 presents the velocity field for the c\=C2 = a = b = \ parameters. Unfortunately, 
there is no traveling wave solution for the temperature distribution T(x,t) = f(() = f(x — 
at 2 ) because the velocity has an explicit time dependence which is transfered to the third 
equation of (1) giving the following expression = —Xf" — 2atf'+v(Qtf which still depends 
on time and cannot be simplified to an ODE of /(C)- 

C. The linear EOS case 

Let's assume now the linear case of EOS where 

p = Ap (37) 

where A is till a constant just to fix the proper physical dimension. From the second case 
of Eq. (1) the velocity field can be obtained for small vs 

v = -^t. (38) 
P 

From the continuity relation the following PDE is available for the density field 

^ - A Pxx = 0. (39) 
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Note, that the former variable transformation t 2 = t 2 /2 is still valid giving us the regular 
Fourier heat conduction or Fick's diffusion equation of 



Pt 2 - Ap x 



0. 



(40) 



The solution is the well known Gaussian function of 



p(x,t) 



-.exp 



—x 



^M~ 2 """V 2At 2 J * V A 6XP V At 2 
According to fl38|) the velocity field is v(x,t) = 2x/t. For the temperature distribution the 
following PDE have to be solved 



-x 



(41) 



X 

Tt + 2— T x = \T X 



(42) 



The solution is very similar to the former case of ( )25|) . Using the same notation the solution 

reads 



, '1 1 3 3 9 \ /l 13 3,. 



For completeness the temperature field is 



T(x,t) = f 



, r , 1 1 3 3x 2 \ x r /l 1 3 3x 2 \ x 
ClM 1 2 " 3 7 ' 2 ' 4At J Vt + 02 U " 3 7 ' 2' 4At J 7!. 



(43) 



(44) 



Figure 3 presents the Kummer M part of the whole solution with c\ = 7 = A = 1, and 
C2 = parameters. The Kummer U part looks very similar. It is straight forward that the 
PDE of the density field (139]) can be investigated with the traveling wave Ansatz fl30l) . too. 
Resulting the next linear ODE of 

Ah" - ah' = 0. (45) 

The solution is obvious h — C\ + c 2 e^ with the density field of p(x,t) — C\ + c 2 e~ a } . 
It turns out from fl38l) that v = —at. After trivial calculations from the heat conduction 
equation we get T(x, t) = c\ + c 2 (x — at 2 / 2). 



D. Investigation of the virial EOS 

We still consider the original system of (JT]), and change the EOS in the second equation. 
Outside the ideal gas and the polytropic equation of state the most important material law 
is the Van der Waals (VdW) formula for gases where p(p, T) = ^ — cp 2 where two of the 
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constants (a, b) can be calculated from the critical point. One of them a is proportional with 
the universal gas constant times the molar mass of the fluid. Plugging this expression back 
to the Euler equation and calculating the derivatives we get constrains to the exponents of 
the Ansatz a, (3, 7. For VdW EOS these constrains are contradictory, so no ODEs can be 
obtained no self-similar solution can be evaluated. 

We may go further and try the virial expansion which can be written in various forms, one 
possibility is the following 

p = ATp(l + Bp + Cp 2 ) (46) 

where A,B,C are constants and can be fixed from experiments. We try the simplest B = 
C = case. Apply the Ansatz of Eq. fl3]). The universality relations are therefore a = 
1, /3 = 7 = 5 — 1/2, which dictate the next ODE system 

~\- \v h ' + 9 h ' + hg' = 0, 

1 1 ,, , Af , Afhl 
— 2 9--^9 +99 = -Af - — 

-f - \vf + 9f' = A/". (47) 

Note, that the first equation is a total derivate and can be integrated immediately which 
fixes g = 77/2 as for the first investigated system in the Eq. (5). So, v(x,t) = ^. (We fixed 
the integration constant to zero.) Inserting this into the second and third equation of ( 14T|) 
we get 

V _ Af , Afh' 
~ 4 ~ h~ : 

-f = Af . (48) 



The solution for / reads: 



f(rj) = cxsin ( -^=j + c 2 cos \^=j ( 49 ) 



which means 



T{x,t) = - t 



x \ fx 
CiSin I — — + C2COS 



<\tj ' \VXt, 

Finally, the formal solution for h{rj) for the density shape function can be derived 



(50) 



-Ac lC os (-fa) + Ac 2 sin + ^ 

KV) = C3 / ■ - , V 7 V dz. (Ol) 

V\A d«n [^ x )+c 2 cos 
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For given parameters like Ci = 0,C2 = c^ = A = X = l the integral can be evaluated in a 
closed form of 



V » 1 



= 3( ' ( ' e :„ +1 (52) 

where i is the complex unit and dilog(x) is the dilogarithm (or Spence's) function defined 
via the following integral 

dilog(x) = — [ — ■ -ds. (53) 



s 

A detailed mathematical description of the properties of the dilogarithm function can be 
found in [251 ] . For completeness we give the final formula of the density 

/ l ,,.'■'/ v ' \ ' 

p(x,t)=r 1 / 2 



(1 _|_ ie^/^Y^il - j e M/V?)^ e -|[-4x/Vt-rf«ofl(l^e fa /^)+diJo ff (l-ie te /^)] 



p2ix/Vi _j_ J 

(54) 

Note, that the velocity and the temperature distributions in the first and the last model are 
very similar to each other. 



III. SUMMARY AND OUTLOOK 



We investigated the basic one-dimension coupled PDE equations describing fluid flow 
with heat conduction and presented self-similar solutions for fluid density, flow velocity and 
temperature. In certain cases we have tried to find solutions for the system of equation 
presented in ([I]). In the case when p ~ p 3 the major achievement is that an exact solution 
is available and the long time decay of density and pressure is a power law of the time 
variable. The other situation when p ~ p 2 is more complicated from the point of view of 
exact solutions, consequently we tried to arrive to results iteratively. In this case solutions 
related to waves has been also obtained with the help of a special function. The following 
case when p ~ p is again a nontrivial situation in the frame of eq. (JT]). At this point one 
can arrive to solutions iteratively, which in certain aspects may resemble to diffusivity in 
dilute systems. Even more for long time behavior a kind of power law decay for density and 
pressure can be also observed. 

Further possible study may be the investigation of the system with other possible 
equations of state, which may be different than the ones considered here. The extension 
of Eq. (JT]) for real fluids(e.g. including the viscous term, temperature or density dependent 
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viscosity etc.) and searching for possible analytic or eventually exact solutions is also a 
possible future problem to be investigated. Regarding computational methods we hope that 
our study can help for benchmark tests for various numerical codes or models. 
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Fig. 1. The density distribution of Eq. 35 for the a = b = ci=C2 = l parameters. 




Fig. 2. The velocity field of Eq. 36 for the a = b = ci = C2 = l parameters. 
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Fig. 3. The temperature field of Eq. 44 only the Kummer M function is presented with 

the 7 = A = Ci = 1 parameters. 



